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Abstract. We prove that any distribution q satisfying the equa- 
tion \7q = divf for some tensor f = (/j), /j £ h r (U) (1 < r < 00) 
-the local Hardy space, q is in h r , and is locally represented by 
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> 
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the sum of singular integrals of fj with Caldcron-Zygmund kernel. 
As a consequence, we prove the existence and the local represen- 
tation of the hydrostatic pressure p (modulo constant) associated 
with incompressible elastic energy-minimizing deformation u satis- 



^ <~| fying |Vu| 2 , |cofVu| 2 £ h 1 . We also derive the system of Euler- 



Lagrange equations for incompressible local minimizers u that are in 
the space if lo ' c (defined in (1.2)); partially resolving a long stand- 
ing problem. For Holder continuous pressure p, we obtain partial 
regularity of area-preserving minimizers. 



1. Introduction 



00 

£<■} . Let O C t", n > 2 be a bounded Lipschitz material body. For Mooney- 

^ ! Rivlin or Neo-Hookean materials [Ba 77], [TO 81], [Og 84], such as vul- 

canized rubber, in the equilibrium state, one is interested in minimizing 
the elastic energy 

(1.1) E[w) := I L{Vw{x))dx, 



n 



for incompressible ll /1,2 -deformations w : Q C M n — > M. n , subject to 
its own boundary condition, and corresponding to a given smooth bulk 
energy L : M nxn —>■ K. Let us define the subspace K 1,r for 1 < r < 00, 
by 

(1.2) i^' r (fi,M n ) := {w G W 1>r (ft,R n ) : cof Vw G L r (fi,M nxn )} , 
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where W 1,r denotes the usual Sobolev spaces (see for example, [GT 97, 
Chapter 7]) and cof P is the cof actor matrix, whose ij-th entry is (— l)*" 1 "- 7 
times the determinant of (n — 1) x (n — 1) submatrix obtained by deleting 
the i-th row and the j'-th column from the n x n matrix P. Using 
the identity P l coiP = Id n det P, it follows that det Vw G L l for any 
w G K 1,2 . Since \P\ = |cof P| for any P G M 2x2 , the function spaces 
K 1,r and W 1,r are equal in M 2 . Let us denote the admissible set of 
deformations 

(1.3) A := {w G J\ 1,2 (n, R n ) : det Vw = 1 a.e. in ft} , 
We call u G A to be a /oca/ minimizer of £?[•] if and only if 

(1.4) E[u] < E[w] for all w G -4 and supp (w — u) C ft . 

Under the hypothesis that the energy density L is smooth, polyconvex 
(convex function of minors) [Ba 77] and satisfies the growth condition 

(1.5) Ci(|X| 2 + | cof X\ 2 ) -C 2 < L(X) < C 3 (l + \X\ 2 + | cof X\ 2 ), 

for all X G M nxn , for some d > 0, C 2 > 0, C 3 > 0, where \X\ 2 := 
trace(X*X), using direct methods in the calculus of variations together 
with weak continuity of the determinant, J. Ball [Ba 77] proved the ex- 
istence of local minimizers u G A of the energy E[-]. An example of 
polyconvex L satisfying the growth condition (1.5) is the stored-energy 
for incompressible isotropic Mooney-Rivlin materials in M. 3 , given by 

(1.6) L(X) = ^(h(X) - 3) + f (J 2 (X) - 3) , 

where h(X) := trace(C) = |X| 2 , I 2 (X) := \ [(trace(C)) 2 - trace(C 2 

| cof X\ 2 , are the first two principle invariants of the right Cauchy- Green 
strain tensor C := X l X and /i l5 fi 2 are positive material constants. 

Though the existence of the local minimizers of E[-\ in A is known for 
over 30 years, the existence of integrable hydrostatic pressure associated 
with such minimizers, the derivation of system of Euler-Lagrange equa- 
tions, and the partial regularity for such minimizers remains a challenging 
open problem. In this article we prove the following results: 

(I) The h r (1 < r < oo) -integrability and local representation 
of any distribution q satisfying the equation Vg = f, where 
f := (/j), fj G h r , the local r -Hardy spaces. (Theorem 2.2) 

(II) The existence of a pressure p G L\ oc if the minimizer is u G 
K^J for some r > 1. (Theorem 3.1) 
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(III) The existence of a pressure p G h 1 if the minimizer u satisfies 
the conditions |Vu| 2 , |cofVu| 2 G h 1 . (Theorem 3.1) 

(IV) The validity of the Euler-Lagrange equations if the minimizer 
is u G K^ c . (Theorem 4.1). The pair (u, p) satisfies the system 

(1.7) div[D£(Vu(a;))-p(a;)cof(Vu(a;))] = in fi, 

where the divergence is taken in each rows. 

(V) The partial regularity of W 1 ' 3 area-preserving minimizers u 
for which the hydrostatic pressure p is Holder continuous with 
exponent < a < 1. (Theorem 5.1) 

The L 2 -version of the result in (I) is classical (see, [Te 01, Remark 
1.4, p 11] ), and plays an important role in incompressible fluids [Te 01]. 
The result in (I) is a crucial ingredient in proving (II) & (III)- The h 1 - 
version of (I) is quite delicate and to the best of our knowledge, it is new 
and may be of independent interest. For the case r > 1, it follows that 
Vq G W~ 1,r , and adapting the classical functional- analytic approach 
demonstrated for r = 2 (see [Te 01], [TO 81]), or arguing directly by 
duality, and solving the equation of the type 

divw = / in V CC U, w = in dV, 

[Ev 98, p. 472-474], one can prove that q G L[ oc . However, both of these 
approaches fail to give informations for the critical case r = 1 and does 
not give a representation of q. Whereas, our unified singular integral 
approach is self-contained, simple and provides the local //-estimate, as 
well as the local representation of q. The main ideas in our proof is 
to represent the localized-mollified distribution q in terms of the New- 
tonian potential in R n and finding its uniform bound in h r , by using 
Calderon-Zygmund estimate [FS 72], [CZ 52]. Finally we show that the 
local representation of q consists the sum of Calderon-Zygmund type 
singular integrals of the tensor f (see equation (2.27) in Section 4). 

For the case n — 2, under the stronger hypothesis that the local mini- 
mizers of E[-] are classical (C^-diffeomorphism), namely in the Sobolev 
space W 2,r for some r > 2, LeTallec and Oden [TO 81] established the 
system of equations in (1.7). For n = 2, Bauman, Owen and Phillips 
[BOP 92] proved that if a minimizer is in W 2,r for some r > 2, then it is 
smooth. For such W 2 ' r , r > 2 minimizers, the authors in [BOP 92] ar- 
gued directly on the level of the Euler-Lagrange equations exploring the 
existence of integrable hydrostatic pressure. Evans and Gariepy [EG 99] 
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proved that any non- degenerate, Lipschitz area-preserving local mini- 
mizers of E[] are C 1,a (f2o), for some < a < 1 for a dense open subset 
Qo C Q. We believe that the Euler-Lagrange equations (1.7) that we 
derived for i^^-minimizers may be useful in understanding the partial 
regularity of such minimizers, as evidenced by the result in (V). 

In order to prove the existence of an integrable pressure p associated 
with the local minimizer u, we only require the additional mild assump- 
tion | Vu| 2 log(2+ 1 Vw| 2 ), |cofVu| 2 log(2+ |cofVu| 2 ) G L\ oc . For n = 2, to 
derive the system of equilibrium equations (1.7) for (u,p) in Q, we need 
u to be in W 1 ' 3 , whereas the best-known previous result in this direction 
were for jy 2 ' r -minimizers for some r > 2. 

We organize the paper as follows. In Section 2 we prove (I); in Section 
3 we prove (II) & (HI); in Section 4 we prove (IV), and finally in Section 

5 we prove (V). Throughout this article C is a generic absolute constant 
depending on n, U, Q, u(f2), 7CC u(f2), r, and L. Its value can vary 
from line to line, but each line is valid with C being a pure positive 
number. 

2. Local integrability of solutions Vg = divf 

We recall some of the basic definitions and terminologies of Hardy 
spaces. Let 1 < r < oo. A distribution / belongs to iJ r (IR n ) if and only 
if / G L r (R n ) and Rj(f) G L r (R n ) (see for example, [St 93, Proposition 
3, p. 123]) for j — 1, • ■ • , n, where Rj is the Riesz transform of / given 
by 

R,(f)(x) := lime f j^f(x - V) dy, c n -- 



■j\j jk^j ■ — -« , , , n+lJ v- yj -y, -« • 7r(n+1)/25 



y\>e W 



so that Rj(f)(£) = iffif- I n short, we will write H r (M. n ) as simply H r . 
For / G H r , the norm is defined as 

\\f\\Hr:=\\f\\ L r+J2mf)hr. 

A standard result [St 70, p. 237] states that a positive function /, the 
Riesz transform Rjf G L\ oc if and only if /log(2 + /) G L\ oc . For 1 < 
r < oo, a classical result asserts that / G H r if and only if / G L r , 
see [St 70, p. 220]. The celebrated Fefferman duality theorem [Fe 71], 
[FS 72, Theorem 2], [St 93, Theorem 1, p. 142] asserts that the dual of H 1 
is the BMO, the functions of bounded mean oscillations. The following 
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theorem is due to Calderon-Zygmond [CZ 52], Stein [St 70, Theorem 3, 
p. 39], and Stein-Fefferman [FS 72, Corollary 1, p. 149-151]. 

Theorem 2.1 (Calderon-Zygmond, Fefferman-Stein). Let 1 < r < 

oo and f G H r . Let G be a C l function on W 1 \ {0} homogeneous of 
degree with mean value over the unit sphere S™^ 1 , that is 



(2.1) / G(x)da(x) = 0. 
Then the function defined as 

(2.2) T f(x) := lim / ^M /( x _ y ) dy 

S ^°J\y\>6 \y\ 

exists a. e. and furthermore, 

(2-3) ||To/||^<CV||/||^. 



In particular, i?j's are bounded linear operator on H r , for any 1 < 
r < oo. Let us recall the definition of local Hardy spaces introduced 
by Goldberg [Go 79]. A distribution / on M. n is said to be in the local 
r-Hardy space, written as / G h r , if and only if the maximal function 

Miocf(x):= sup \{p £ *f){x)\ 

0<e<l 

is in V , where p e := e~ n p(x/e), is a standard approximation of the 
identity. The h r norm of / is defined to be the L r norm of the maximal 
function A4\ oc f. It follows that if / G h r then rjf G h r for any smooth 
cut-off function and H r C h r . For bounded Lipschitz domain Q C M n , 
we adopt the definition of Hardy spaces h r (Q) introduced by Miyachi 
[Mi 90]. A distribution / on Q is said to be in h r (Q) if / is the restriction 
to fi of a distribution F in ^(M 71 ), i.e., 

h r (Q) ■= {/ g X>'(f2) : 3F G /i r (R n ), such that F| n = /} 
= h r (WL n )/{F G /i r (M n ) : F = on fi}. 

The norm on this space is the quotient norm: the infimum of h r norms 
of all possible extensions of / in M n . For 1 < r < oo the spaces h r (Q) is 
equivalent to L r (Q). For smooth bounded domains Q, the Theorem 2.1 
is valid for / G h l (n), see [Mi 90], [CKS 93]. 
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Theorem 2.2. Let U C W 1 , n > 2 be a bounded Lipschitz domain and 
1 < r < oo. Let f = (/j) suc/i i/iai /j G h r {U), for 1 < z, j < n. T/ien 
£/te distribution q : C^°(U) — > K. defined by 

(2.4) Vg = div/ ^^ (Vg, v) = - f{x) : Vv(i) dx 

/or v G Co°(C/, M n ), is m /i r (^), /or any V CC C/ w/iere A : B :— 
trace(A t B) = Y^^albp for A, B G M nxn . Furthermore, q is locally 
represented by sum of singular integrals of /j (see equation (2.27)), and 
for any V CC U, there exists C > 0, depending only on U, V and r such 
that 

\\q\\hr(V) < C||/IU r (v) • 



Proof of Theorem 2.2. Let U C M n , n > 2 be a Lipschitz domain. Let 

f : = (/j) g M nxn and /j G /i r (£7), for 1 < r < oo and 1 < i, j < n. Let 
g G £>'(£/), such that 

(2.5) Vg = divf in V'(U) . 

Our idea is to mollify the equations in (2.5) and obtain uniform bound for 
the mollified q, by using Calderon-Zygmund estimate. Let V CC U be 
a sub-domain and < e < dist(V, dU). Let p £ be the usual mollification 
kernel, and define convolution q £ : V — > R by 

q £ (x) = (q*p £ )(x) := (q, (p e ) x ) for x G V, where {p e )x{y) ■= p £ {y-x), y EU 

Then by the standard properties of the mollification [DL 88, Proposition 
1, p492], q £ is smooth and for any 1 < % < n 

d _ dq_^ __ dp £ 

Hence mollifying the system of equations in (2.5), we obtain 

(2.6) Vg e = divf £ in V, 



where the divergence is taken in each rows of f e := ( (fj) £ I , and (fj) £ '■= 
fj * p £ is the mollification of f. Since /j G h r (U), we conclude that 

(2.7) (/j) e -/j strongly in h r (V) as e - 0, 

for all 1 < i,j < n. Applying the divergence operator to the both sides 
of the above equation, we obtain 

(2.8) Ag e = div(divf e ) in V. 
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Since there is no control on the boundary values, we need to localize the 
equation (2.8). Let W CC V CC U. Let rj G C °°(M n ), < 77 < 1 be a 
cut-off function such that r] = l'm.W and 77 = outside V. Let g £ := rjq £ 
be the localization of q £ . Then g e is the solution of Poisson equation 

(2.9) Aq £ = f £ in R n , 
where 

(2.10) f e := r/Ag e + 2(Vg £ , Vr?) + q £ Ar, 

= r\ div(div f e ) + 2(div f £ , Vr?) + g £ A77. 

Therefore q £ is represented by the Newtonian potential of in W 1 . In other 
words, 

(2.11) q e (x) = ~ I Hx-y)f E (y)dy, 

where $ is fundamental solution of the Laplace equation in R n and is 
given by 

{ — ^: log \x\ if n = 2 

n(n-2)a(n) M™" 2 " U — * ' 

n/2 

for x G IR n \ {0}, and a(n) := r r"n +1 -, is the volume of the unit ball in M. n . 
Using (2.10) in (2.11), we obtain 

(2.13) q £ (x) = - ri(y)$(x-y)dw(divQdy 

Jw. n 

+2/ ((div f £ ,X7ri) + q £ Ari)$(x-y)dy 
Jm. n 

:= -l}(x)-2I?(x)-I?(x), 



r)(y) $(x - y) div(divf e (j/)) dy 

(divf e (y),V?7(y)) $(x-y)dy 
Qe(y) &ri(y) $(x - 1/) dy 



where 




(2.14) 


£(*) 


(2.15) 


4 2 (*) 


(2.16) 


Ilix) 
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By direct computations, observe that, for 1 < i, j < n 

1 VVi 



(2-17) (77$) = rjyMy) 



u n \y\ n 



(2.18) („*V = %w 4(y) 



3 uj n \y\ 

1 fx ViVi\ V 

1 da - n 



u n v \y\ J \y\ 

where 5ij is the Kronecker delta and u n := na n is the surface area of 
the unit sphere § n_1 . We now establish an uniform local /i r -estimates 
(1 < r < 00) for q e through the following steps. 

Step 1: Limit of i|. Let us fix x G W CC V CC 17. Since A77 = 
on VK, the integrand in 7g(x) is smooth. Since q e is determined up to 

a constant, we can add a constant to y h-> A?7(y)$|x — j/|, if nessecary, 
to ensure that it has vanishing integral. For each fixed x G W, let 
v x : V^ — > W 1 be the solution of the Dirichlet problem 

(2 19) { divv ' x ^ = &r}(y)$(x - y) for y G V 

^ ' ' I Va; = on dV . 

Then using (2.19), integrating by parts, and the convergence of f £ in 
(2.16), we obtain 

(2.20) I 3 £ (x) = [ q £ (y)Ar ] (y)$(x-y)dy 

q E (y)divv x (y)dy 

- I (Vqe(y),v x (y))dx 

(div f £ (y),v x (y))dy 

f e (y) : V y v x (y)) dy 

1 

%) : V y v x (y) dy as e -> 

:= Iq(x) ioi x eW CC V. 

Since f e — ► f strongly in /i r (V,M nxri ), it follows that J^ — > Jq strongly in 
/i r (iy). 



EULER-LAGRANGE EQUATIONS 9 

Step 2: Limit of I 2 . Let us fix x G W CC V CC U. Integrating by 
parts, invoking (2.17) and letting e->0we have 

(2.21) I 2 (x) = f (divUy),Hx-y)Vv(y))dy 
= - [ f £ :V y U{x-y)V7])dy 



ft: U{x 

JR n V 


- y) V 2 r/ - 


(y — x) <S> V?? N 
uj n \y-x\ n , 


) dy 


- / f : ($(x- 

JR n V 


- y) V 2 r/ - 


(y — x) ® V??^ 
u n \y-x\ n ) 


\dy 


ll(x) xeW. 









Using the strong convergence of f e in h r (V), again it follows that I 2 — > Iq 
in h r (W). 

Step 3: Limit of I]. Integrating by parts twice the integral in (2.14) 
and using (2.18) 



l]{x) = / divdivf £ (y) r)(y)$(x - y) dy 

JR n 

= [ Uy)--V 2 y (v(yMx- y ))dy 



Uv) : fax - y) V*v) - ± ^»(»-') + (»-»)«^ ) dv 

V w n \x-y\ n } 

-- f «*,) : (li. - » (»-*) «<»-*) ) _L) rfy 

^n y R » V \x-y\ J \x-y\ J 

:=ll 1 {x) + ll 2 {x), XEW, 

where Id n is the n x n identity matrix. Using the convergence of f £ , 
observe that as e — > 0, 

riif„\. f ? . (*(- „\V72„ Vrt®(y-x) + (y-x)®Vr] 



I\x):= f £ :(Hx-y)V 2 V -^^ V , w ' ' ) dy 

V u n \x-y\ n J 

f : U{x - y) v 2 v - vq®(y-*) + (y-*)®V'A d 

(2.22) :=/ n (x), s£W. 

In order to estimate I] 2 , define the kernels fiy : lR n \ {0} — > K by 

(2.23) %(</):=c^-n^, |/el B \{0}, i,J = 1, 
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Since na n = u n , integrating by parts, observe that for any i, j = 1, • • • , n, 



Slij(y) da(y) = / ((% - ny^j) da(y) 



[Or, 



nSij - n / y { yj da(y) 



n-l 







uj n 5ij -n T^—yi dy 
Jb 1 dy 3 

u n 5ij - n5ija n 
0. 



Hence each f^- satisfies all the conditions of Calderon-Zygmund Kernel 
[St 70]. Therefore, 

(2.24) /»(,) : = -J- / ,f. : U - „ <» = ?® %Zf> ) » ) 

is the sum of Calderon-Zygmund type singular integrals with the homo- 
geneous kernel fiy. Since f G h r (U,M. nxn ), 1 < r < oo, by Theorem 2.1 
J 12 G /i r (jy). Furthermore, the following sum of singular integrals 

1 f ( (y - x) ® (y - x) \ dy 



^Jm« V \x-yr J 



\x — y\ n 



exists for almost every x G W CC V and is in h r (W). From the singular 
integrals (2.24) and (2.25), by Theorem 2.1, we have 



I?(x)-I?(x) 



U n . ■_-. JM n \ J \ X ~ y\ 

Hence there exists C := C(V, W,r) > such that 

(2.26) ||/ £ 12 - Il 2 \\ h r {w) < C*]T \\U])e ~ f]\Wiy) ^0 as e - 0. 

3=1 

Step 4: Explicit representation of q. To complete the proof, let us 
define the potential q : W — > M by 

g(x) := -(/ n (x) + J 12 (x) + 2J 2 (x) + / 3 (x)) . 
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Then from (2.20), (2.21), (2.22), and (2.26), we conclude that q £ -+ q 
strongly in h\ oc for any 1 < r < oo, and hence and q is represented as 



(2.27) q(x) = f f : ($(x - y) V 2 V - V y v x ) dy 
Ju 



dy 



x-y\ 
dy 



■ / f : V?7 ® (y — x) — (y — x) ® Vr/ 
+ — r)f:[Id n - n 

u n Ju V |x-yr / \ x -y\ 

for any a; G W. Since g is the strong limit of the family q e in W, it is 
independent of the choice of the cut-off function 7]. This completes the 
proof of Theorem 1.1. □ 



3. First Variation of Energy and the existence of 
hydrostatic pressure 

Let O C M", n > 2 be a smooth, simply connected and bounded 
domain and let L : M nxn — > M. be smooth function. We are now in 
a position to establish the existence of integrable hydrostatic pressure 
associated with the local minimizers of the energy 



E[w] := / L(Vw(x))dx. 



(3.1) 



for incompressible jy 1,2 -deformations w : Q C lR n — > M. n . By direct 
computation, observe that Mooney-Rivlin bulk-energy given by 



(3.2) 

satisfies the following. 



L(X) = ^(|Vu| 2 -3) + ^(|cofVu| 2 -3) 



DL = HiP + fj, 2 



( coi{SQ)\ : {SP)\ -coi{SQ)\ : {SQ)\ coi{SQ)\ : {SP)\ \ 
-coi{SQ)\ : {SP)\ coi{SQf 2 : {SP)\ -coi{SQ)l : (SP)\ 



\ coi(SQ)l : (SP)l -coi(SQf 2 : (SP) 3 2 coi(SQ)l : (SP)* J 

where Q := cofP, and (SX)j is the 2x2 submatrix obtained by deleting 
the i-th row and the j-th column of the matrix X G M 3x3 . Furthermore, 
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the Cauchy- Green strain tensor is given by 

(\Q2\ 2 + \Qs\ 2 -<Qi,Q 2 > -(Qi,Q s ) \ 

-<Qi,Q 2 > \Qi\ 2 + \Q 3 \ 2 -<Q 2 ,Q 3 > 
V -<Qi,Q 2 > -<Q 2 ,g 3 > \Qi\ 2 + \Q2\ 2 J 



(DL(P)) t P = fx 1 P t P + fi 2 



for all P G M 3x3 , where Q t := (cofP); := ((cofP)j, (cofP) l 2 , (cofP)|) be 
the i-th row of cofP, i = 1,2,3. Motivated by the above calculations, 
assume that L satisfies the following growth condition. 

(3.3) max (\L(P)\, \{DL{P)) t P\ J < C(l + |P| 2 + |cofP| 2 ), 

for some C > 0, for any P e M nxn . 

Now we prove the existence of an integrable hydrostatic pressure q on 
the deformed domain u(Q) and establish an explicit representation of 
the pressure q in terms of Calderon-Zygmund type singular integrals of 
the Cauchy-Green strain a := (PL(Vu))*Vu) o u _1 in u(f2). Our proof 
consists of deriving the first variation of the energy E[-], obtaining the 
equation Vg = diva, and then finally use Theorem 2.2. 

Theorem 3.1. Let L : M nxn — > K be smooth and satisfies the growth 
condition (3.3). Assume that we A be a continuous and injective local 
minimizer of E[-], such that |Vtt| 2 , |cofVu| 2 e hl oc (Q) for some 1 < r < 
oo. Then there exists a scalar function q G h"[ oc (u(Q)), such that 

\\q\\ h r ( v) < C (\\\Vu\ 2 \\ hr{u _ 1{v)) + |||cofVw| 2 ||, r( ^ 1(y)) ) , V CC u(Q), 

for some C > (depending on r, V, n and u{Vt)) and the pair (u,q) 
satisfies the integral identity 

(3.4) / DL(Vu(x)) : V(vou)dx= / q(y) div v(y) dy 

Jq Ju(n) 

for all v E C^°(u(n),R n ), where A : B := tr(A*B) = £y=i ajbj is the 
scalar product on M nxn . 
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Remark 3.2. Let W CC V CC u(fi), and 77 G C£°(F) be a cut-off 
function such that rj = 1 on H 7 . Then g is locally represented as 

(3.5) g(x) = f a : ($(x - y)V 2 r] - V y v x ) dy 
Jv 

H / a : ( V?? (g> (y — 2) — (y — x) <S> Vr/ I — 

unJv V /f-j/P 

1 f _ ( (y - x) ® (y - x)\ dy 



^n Jy V \x-y\ 2 J \x-y\ n 

for any x G W, where $ is Newtonian potential in M. n defined in (2.12) 
and v x as defined in (2.19). 

Remark 3.3. In the study of regularity of finite energy deformations, 
Sverak [Sv 88] proved that for any W 1,n -deformation w with det Vw(x) > 
0, a.e., there exists a continuous function uj on K with uj(0) = such that 

|w(x) — w(y)| < cu(|x — y\), for any i,j/6S] CC M n . 

It is also well-known any W^ 1,n -deformation w for which the distortion 
function K(-,w) := | Vw(-)| n /det Vw(-) G U for some r > n — 1, is 
a homeomorphism. Thus in particular, area-preserving VK 1,r (r > 2)- 
deformations in the plane are continuous and open maps. However, in 
general for n > 3, any deformation w G K ,2 may be totally discontinu- 
ous, see [Sv 88, p. 119]. 

In order to prove Theorem 3.1, we establish the following first variation 
of the energy integral E[-}. 

Lemma 3.4. First Variation. Let uE Abe a local minimizer of E[-}. 
We further assume that u is a continuous and an injective map. Then u 
satisfies the following integral identity 

(3.6) DL(Vu{x)):V{vou)(x)dx = 0, 
Jn 

for all smooth, compactly supported and divergence free vector fields v on 
u(Q). 

Proof: By the invariance of domain u(Q) is open and u : Q — > u(Q) is a 
homeomorphism. Let v G C^°(u(f2), MJ 1 ) be a vector field with div v = 0. 
For each y G u(f2), consider the unique smooth flow (p(y, •) : 1R — * u(f2) 
given by 

(3.7) ^(l/,*)=v(0(i/,t)) in R, <j>(y,0) = y. 
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d 
Using the relations -^-^detP = (cofP)* and P (cof P) 4 = Id n detP, by 

3 

a direct calculations we observe that 



(3.i 



d 

dt 



(det V y (/)(y, t)) = det V y <f>(y, t) div v = 0. 



Since det V y (j)(y, 0) = 1, from (3.8) it follows that det V y (j)(y, t) = 1 for 
all £ G R and y £ u(Q). Consider the map w : f2 x R — ► u(f2) defined by 

w(x, £) := 0(-, t) o u (x) = 0(u(x), t) for any £ e R, iffi. 

Let y := suppv C u(fi), then v(u(x)) = for u(x) ^ V. This in 
conjunction with the uniqueness of implies that 0(u(x), £) = u(x) for all 
points x such that u(x) $_ V . Since Q is bounded, u is continuous and V 
is compact, £)' = u~ 1 (V) is a compact subset of Q. Hence supp(w(x, t) — 
u(x)) C Q'. Furthermore, detV x w(x, t) = det V y (j)(y, t) det Vu(x) = 1. 
Therefore, w(-,t) e „4 and supp(u — w(-,£)) C Q for all <6l. Since u 
is a local minimizer of E[-], 

E[u] < P[w(-,t)] for all t e R. 

Thus in particular, 

0= — / L(Vw(x,t))dx 



i=0 



E 

2 

E 

M=l 
2 

E 

2 

E 



L}(Vw(x,£))^f|^-(x,t))dx 



t=0 



5 



t}(Vw(x,«))^^(» 



i(x),t) I dx 



t=o 



5 



L}(Vw(x,£)) — - (^(0(w(x),£)) dx 



dxi 



L}(Vu(x)) A (</(u(x))) dx 



i=0 



/ DL(Vu(x)) : V(vou)(i)di, 
in 



for all smooth, compactly supported and divergence free vector fields on 

dL 

u(f2), where L)(P) := t-^(P) • This proves the Theorem. □ 

J op l j 
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Proof of Theorem 3.1: Let 1 < r < oo and U' CC U. Let u E A 
be a local minimizer of E[-] such that |Vu| 2 G h r and |cof Vu| 2 G h r (U') 
for some 1 < r < oo. Assume further that u : Q — > u(Q) is continuous 
and bijective map. 

Now define g = (g 1 , ■ ■ ■ , g n ) : C7 1 (u(fi), R n ) -> R by 
(3.9) (g, v) := / DI(Vu(i)) : V(v o u)(x) dx, 



for all v = (v 1 , ■ • • , v n ) G Qj (u(fi), W 1 ). In view of the volume constraint 
and growth condition (3.3), it follows that 

(3.10) |(g, v)| < C{\ + ||Vu|| L 2 (n) + ||cofVn|| L2(f7) ) ||Vv|| LO o (u(n)) , 

for any v G Cg(u(f2),IR n ). Hence g is a continuous linear functional on 
Cg(u(f2),IR n ). Using the the first variation (3.6), we conclude that 

(3.11) (g,v) = Vve^(u(0),r), divv = 0. 

Hence there exists q G V'{\i{Q)) ( see [Te 01, Proposition 1.1, plO]), such 
that 

(3.12) g = -Vg inP'(u(n),R n ) 

modulo translation of a constant. In order to obtain h r estimates of q, 
for 1 < i,j < n, let us define a^ : £1 — > R by 

(3.13) <r){x) := Y, L U Vu ( x )) «^( x ) for x e fi ' 

fc=i Xfc 

so that, the Cauchy- Green strain tensor on Q is given by 

(3.14) a := (a}) = (£>L(Vu))'Vu 

Define the r/'-th component of the Cauchy- Green Strain tensor a!- on the 
deformed domain u(Q) by 

(3.15) 5-]- := (Tj o u _1 on u(fl), i, j = 1, • • • ,ra. 

The growth condition [cr*-| < C(|Vu| 2 + |cofVu| 2 ) and |Vu| 2 , |cofVu| 2 G 
LlogL yields a) G /^(K). If u G K^ r {n,R n ), 1 < r < oo, from the 
definition of <r!-, a*-, and the condition (3.3) on L, it follows that 



(3.i6) / icsjr = / k}r 



C C IIVu||/ 2 r( u -l(y)) + ||cofVu|| L , 2r( u -l(y)) 
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for any V CC u(f2). Therefore, if |Vu| 2 G h r and |cof Vu| 2 G h\ oc for 
some 1 < r < oo, from (3.16), we have 

a ■= (a)) G h r loc (n,M nxn ) and a : = (crj) G /i[ oc (u(fi), M nxn ). 

Observe that, from the definition of g in (3.9), a!- in (3.13), a* in (3.15), 
and change of variables, 

(3.17) <g,v)= J2 [ Ll(Vu(x))^-^ou)(x)dx 

= J2 Ll(Vu(x))— (u(x))^-(x) dx 
i,j,k=i J ° dy i dXk 

<j{x) : V u v(u(x))rfx 
a(y) : Vv(y) dy 

r u(0) 

= — ( diver, v) 
for any v G C 1 (u(fi),M n ). Hence 

(3.18) g = -diva in£>'(u(fi),M nxn ) 

where the divergence is taken in each rows. Therefore, combining (3.12) 
and (3.18), we get 

(3.19) Vg = diva in V'{u{tt),M nxn ). 

By taking f = a, and U — V CC u(Q) in (3.19), from Theorem 2.2, we 
conclude that q G /?{ oc (u(f2)), it satisfies the local representation (3.5), 
and 

(3.20) \\q\\ h r {v) < C\\a\\ hr{v) 

< C (lllVul 2 !!^-!^)) + |||cofVu| 2 || h r (u ~i (K)) ) , 

for any V CC u(f2), for some C > 0, depending on r, V, n and u(f2). 
Since q G L\ oc , from (3.12), it follows that 

(g,v) = -(Vg,v) = (g,divv) = / q(y) divv(y) dy. 

iu(!i) 
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for any v E Cg(u(fi), R n ). Hence 

(3.21) / DL(Vu(x)) : V(vou)(i)di= / g(y) div v(y) dy, 

JU Ju{Q) 

for any v E Cg(u(ft),IR n ). This proves the Theorem. □ 

4. Derivation of Euler-Lagrange Equations 

Theorem 4.1. Lei ft C K n , n > 2, &e a smooth, simply connected 
and bounded domain. Let u E An -ftT lo '*(ft, R n ) for some s > 3 be a 
continuous and injective local minimizer of E[-}. Then the hydrostatic 
pressure p := q o u E L^ c (ft), and the pair (u,p) satisfies 

(4.1) / DL(Vu(x)) : V(j)(x)dx= / p(x) coi (V u(x)) : V0(x)cfe, 
Jn Jn 

for all (f) E CQ(f2,IR n ) ; where q E L^ c (li(ft)) as in Theorem 3.1. In other 
words, the pair (u,p) satisfies the system of Euler-Lagrange equations 

div [DL(Vu(x)) - p(x) cof (Vw(x))] =0 in ft, 

in the sense of distribution, where the divergence is taken in each rows. 

Proof. Let O C 1" be a smooth, simply connected domain. Recall 
that K 1 ' 3 := {w E W l > s : cof Vw E L s } and A := {w E K^ 2 (Q,R n ) : 
det Vw = 1 a.e.}. Let u E A D if lo '*(ft, R n ), s > 3 be a continuous 
injective local minimizer of the functional E[-}. By Theorem 3.1, there 

exists q E L^ c such that the pair (u, q) satisfies the identity (3.21). Let 
u _1 : u(ft) — > ft be the inverse of u. Then using the volume-constraint 
we obtain 

VyU^iy) = (V x u(x)) _1 = (cofVu(x))*, y = u(x), 

and hence by the change of variables 

/ \Vu~ 1 (y)\ 2 dy = / |cof Vu{x)\ 2 dx < oo. 
Ju(n) Jn 

Using the relation coi(XY) = cof X cof Y, for X, Y E MP 1 * 71 , observe 
that 

Id n = cof (VyU^ 1 Vu) = cof (V^u" 1 ) cof (Vu) = cof (V^lT 1 ) (Vu) _t , 

and hence 

co^VlT 1 ) =(Vuf. 

Since u E if£'(fi,R n ), it follows that u" 1 E K^ s c (u(Q), ft) for s > 3. Let 
V CC u(ft) and € C^(u _1 (F),R re ). Then the composition o lT 1 E 
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Wo' s (V,R n ). Hence there exists v e G C%(V,R n ) such that v £ -> ^ : = 
0OU- 1 strongly in W M (V;M n ) as £ -» 0. Let [/ := u _1 (y). Then Holder 
inequality yields 

DL(Vu) : (v(v £ ou) - V(tpou)jdx 

/ "(Vu)'DL(Vu) : ( V,v e (u) - V 2 ^(u) jdx 



< C|I Vu IIl2"'(C/) l|V(v e -y)\\L*(y), 

where s' := s/(s — 1). Notice that s > 3 yields 2s' < s and hence 
Vu G Lf oc (tt) C Lj^fi). Therefore, from (3.9) we obtain 



(4.2) (g, v £ ) = / DL(Vu(x)) : V(v e o u)(x) dx 

/ DL(Vu(x)) : V(0 o u 1 ou)(x)dx ase^O 

DL(Vu(x)) : V0(x)rfx. 



(V) 

(V) 

u-i(V) 

Since Vu, cof Vu G Lf oc , q G Lj^ and L^ c 2 C L\^~ for s > 3, applying 
change of variables in (3.21), and letting e — > we obtain 

(4.3) (g,v e ) = / q(y) trace (Vv e (y)) dy. 

g(u(x)) trace! V u v £ (u(x)) J dy 
/ q(u(x)) trace! V(v e o u)(x) (cof Vu(i))' 1 dx 

Ju-i(V) V / 

q(u(x)) cof (Vu(i)) : V(v e o u)(x) dx, 
/ q(u(x)) cof(Vu(x)) : V(0 o u~ 1 ou)(x)dx 

iu-!(V) 

q(u(x)) cof (Vu(x)) : V</>(x) dx . 

'u-i(V) 

Hence from (4.2) and (4.3) we obtain 
/ DL(Vu(x)) : V0(x) dx= <?(u(x)) cof (Vu(x)) : V0(x) dx , 
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for any G Co(u -1 (V), M"). Finally choose a sequence of smooth, simply 
connected sets 14 CC 14+i CC u(fi) sub-domains such that u(fi) = 
UfLiVk- Utilizing the foregoing arguments, there exists qu E L s / 2 (14), 
k > 1 such that 



(4.4) / DL(Vu) : V0 = / g fc (u) cof (Vu) : V0, 

for G C*Q(u~ 1 (l4),lR n ). Since u is locally volume-preserving homeo- 
morphism, fi = U^ 1 u~ 1 (Vfc) is an open covering of fi and u _1 (14) CC 
u _1 (14 + i). Using the identity div cof Vu(x) = and invertibility of 
Vu(i), from (4.4) it follows that q^ is unique up to a translation of a 
constant. Thus adding constant terms as necessary to each qk, we deduce 
from (4.4) that for each fixed k > 1 

qi(z) = qk{z) for z E V h \<i< k. 

We finally define q : u(fi) — > K. as g(z) := g/c(z), for z G 14, so that 
g G L^ c (u(fl)). This proves that for any G CQ(fi,IR n ), the pair (u, g) 
satisfies 



DL(Vu(x)) : V0(x) dx = / q(u(x)) cof (Vu(x)) : V0(x) dx . 

Now let us define the pressure p on fi by 

p(x) := q(u(x)) for x G fi. 
Then for any fc > 1, 

\p(x)\ s/2 = [ \q(u(x))\ s/2 dx = [ \q(z)\ s/2 dz < oo, 



and hence p E L^ c (fi) and the pair (u, p) satisfies 

(4.5) / DL(Vu(x)) : V0(x) dx = p(x) cof (Vu(i)) : V<j>(x) dx , 
./n Jn 

for any G C*o(fi,IR n ). In other words, (u,p) satisfies the system of 
Euler-Lagrange equations 

div [DL(Vu{x)) - p(x) cof (Vu(x))] = 0, in fi. 

in the sense of (4.5). This completes the proof. □ 
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5. Partial Regularity of area-preserving minimizers 

For n = 2, as a consequence of the Euler-Lagrange equations (1.7), 
together with the standard elliptic estimates [GM 79], we establish the 
following theorem. 

Theorem 5.1. Let Q C K 2 be a smooth, bounded simply connected do- 
main and let L : M 2x2 — >• M be smooth, uniformly convex, such that DL 
has linear growth and D 2 L is bounded. Let u G W 1,3 (Q,M. 2 ) be an area- 
preserving minimizer of the energy E[-]. Furthermore, assume that the 
associated hydrostatic pressure q on the deformed domain u(Q) is C 0,a 
for some < a < 1. Then Vu are Holder continuous on a dense open 
set Qo C Q. 

Proof. Since u G W 1,3 (Q,M. 2 ) and u is area-preserving, u(f2) is open 
and u is a homeomorphism from Q to u(f2). By Theorem 4.1, there exists 
q G P loc (u(f2)) and the pair (u, gou) satisfies the system 

(s- 1 ) E ^- (||( Vu ) - p ^ ( cof Vu )j) = °' in ^ l = *> 2 ' 

where p := q o u. Assume that q G C°' a (u(Q)). Since u G W 1,3 , Sobolev 
imbedding theorem yields u G C 1//3 , and hence p(x) = q(u(x)) is Holder 
continuous with the exponent a/3. Let F : Q x M 2x2 — > R be the 
free-energy defined as 



[ 2x2 , 



F(x,P):=L(P)-p(x)detP xett, Pe 
so that we can rewrite the nonlinear system (5.1) as 

2 r, 

( 5 - 2 ) E^( A 5^ Vu )) = ' in ^ z = 1 ' 2 ' 

where 

f)F f)T 

*faP) ■= W ^ p ) = W iP)-v{x){coiP)). 

Lett/ CC n. Since |cofP| = |P| for any P G M 2x2 , |DL(P)| < C(1 + |P|) 
and D 2 L(P) is bounded, 



(5.3) \A)(x,P)\<C(l + \P\), 



dA\ 



<c, 
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for any x G U, P G M 2x2 . By Holder continuity of p, it follows that 

\A*,(x,P)-AUy.P)\ l(cofP)!-| 

(5-4) ' A \ - ln r ' l = \p(x) - p(y) | ' 



P\ " v ' rya " l + \P\ 

<C\x-yf\ 

for any x G U, P G M 2x2 . By direct calculations and the ellipticity of L 
it follows that 



dA) d 2 F 

—^(X,P)^kl --■■ g—, 

d 2 L 



( 5 - 5 ) idr(*> p )^ = i+rfaPfafa 



k (P)^ki-2p(x)det£ 

OPjPf 



|2 



> A K| 2 -2p(x)dete 

:= 1(3,0, for P=( P ;),e=(6 



rr2x2 



where Ao > is the ellipticity constant of L. Completing squares, observe 
that 

(5.6) £M)=|^-2^deU 

^0 <^0 

— £ll + £l2 + §21 + ^22 — 2— (^11^22 — §12§2l) 

61 _ ~j~& 2 j + ( ^12 — T~§21 



Ao 



Similarly, we 


obtain 












(5.7) 


Ao 


= (e 22 - ^61 


)\(u 


-i» 


"♦(- 


p 2 \ 
A^ 


(e? 1 +e 1 2 : 


Addii 


ig the identities (5.6) and (5.7), 1 


we obtain 








(5.8) 


I 

2 ~r 


( 6i " i^y 


+ U12 ~ 


\ 2 












+ 1 62 " i^\ 


) 2+ («- 


-£<■) 


"♦(■ 


P 2 ~ 


)iei 2 




> 


kk 
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Thus from (5.6) and (5.8), it follows that the map P *—> A(-, P) is strongly 
elliptic if there exists [Aq > such that 

_*(x,P)£^ >-±h-Z-\ |e| 2 > /i |e| 2 , for xen,P,(6 M 2x2 , 

which is equivalent to assume that 

(5.9) P 2 < A 2 - 2A /io =► (p - /io) 2 < (A - Po) 2 - 

Since p is defined up to addition of arbitrary constant, thus the inequality 

(5.9) is satisfied in subdomain U CC Q if and only if 

(5.10) oscup < A . 

Since p is Holder continuous, the estimate (5.10) holds for any subdomain 
f/cO with sufficiently small diameter. Hence A(x, P) is strongly elliptic 
in P for each i e [/ CC fi, for sufficiently small diameter. This proves 
that A 1 , (x, P) satisfies all the conditions of Giaquinta-Modica in [GM 79] 
on [/ CC H, with diameter of U being small. Hence by [GM 79, Theorem 
1], we conclude that Vu is Holder continuous on a dense open subset Uq 
of U. By standard covering arguments we conclude the proof. □ 
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